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The Schwinger-Dyson, Bethe-Salpeter system of equations are the link between coloured quarks and gluons, 
and colourless hadrons and their properties. This talk reviews some aspects of these studies from the infrared 
behaviour of ghosts to the prediction of electromagnetic form-factors. 
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I. FERMIONS, BOSONS IN THE CONTINUUM 

Strong coupling QCD in the continuum is the physics of 
the real world of gluons, ghosts and quarks: a world in which 
quarks have their physical masses and the pion at 140 MeV 
is by far the lightest hadron. This is the world that can be 
studied using the field equations of the theory and the con- 
sequent bound state equations: the Schwinger-Dyson and 
Bethe-Salpeter equations iQl-dl]. These constitute an infinite 
set of nested integral equations that we cannot solve except in 
rather particular truncations. The best known approximation 
is of course that of perturbation theory, where every Green' s 
function is expanded in powers of the coupling and we stop 
at some finite order. However, the key physics that builds 
hadrons and determines their properties is one of strong cou- 
pling. Nevertheless, the beauty of the perturbative expansion 
of a gauge theory, whether QED or QCD, is that the Green's 
functions are multiplicatively renormalizable and fulfil the re- 
quirements of gauge invariance at every order of truncation. 
Indeed, these are key properties of the full theory for any 
strength of coupling. They thus serve as guides for how we 
might usefully truncate a strong coupling expansion. 

To start the strong coupling study of bound states one 
needs to know their basic building blocks, the propagators for 
quarks, gluons and ghosts and of course their interactions. The 
fact that the field equations are an infinite set of coupled equa- 
tions may make one think that even to solve for the propaga- 
tors (the two-point functions) one needs to know not just the 
three and four point interactions, but even the 24-point func- 
tion and beyond. If that were the case we could hardly make 
progress. While the fermion and boson propagators depend on 
two functions that multiply the two independent spinors and 
tensors, respectively, the fermion-boson interaction depends 
on 12 functions and the 10-point function some huge number 
This looks impossible. Fortunately, the two point functions do 
not depend on each of these independently, but only collec- 
tively. The properties of multiplicative renormalizability and 
gauge invariance of the two point functions pull through to 
these Green' s functions the key elements of the higher point 
functions. 

In general, the inverse boson propagator carrying momen- 
tum p, as shown in Fig. 1 for QED, has two independent tensor 
structures, so 



q=k- p 



q= k-p 



p 



FIG. 1: Schwinger-Dyson equations for the inverse boson and 
fermion propagators in QED. The dots denoted fully dressed quanti- 
ties 



However in a covariant gauge the renormaUzation functions 
A and B are not independent. Indeed, in the Landau gauge 
A(p^) - B{p^). For the boson propagator in QED, two con- 
ditions must be imposed to ensure this. The first is that the 
full fermion-boson vertex satisfies the Ward-Green-Takahashi 
(WGT) identity. To see how this works [4] consider the 
Schwinger-Dyson equation for the inverse photon propagator. 
Fig. 1: 



n^^\p) = Ifjip) 

I— r d'k Tr [/ S f(k) r(k, q) S F(q)] 



(2) 



(27, 



where q = k - p . We contract the propagator with p^py, and 
implementing the WGT identity 



q.r'^dcp) = Sp\k) - Sp'ip) 



yields 



p.n'-'py = [ d'k Tr [p (Sfik - p) -Spik))] 

{2.7tr J 



(3) 



(4) 



IV^\p) = A{p^)p^g''^ - B{p^){fp^ 



(1) 



since the bare propagator is transverse in the Landau gauge. 
The loop integral in Eq. (4) runs over all momentum com- 
ponents from -oo to +00. If the integral were convergent, it 
would be obvious shifting the first term so k - p ^ k that on 
integration the answer would be zero. However, with a cut-off 
regulator this is no longer the case. We have to ensure that the 
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integrals are regulated in a translationally invariant way if the 
integral is to vanish. Thus one must not just satisfy the WGT 
identity, but to be able to make the infinite integrals finite im- 
pose multiplicative renormalizablility too. 

It is in fact only very recently that we have learnt how to 
construct a full fermion-boson vertex that ensures both the 
fermion and boson propagators in QED are multiplicatively 
renormalizable |5]. An illustration of this is shown in Fig. 2 
for the photon renormalization function [6] . (The boson renor- 
malization functions plotted here, and later for QCD, are all 
defined as the propagator multiplied by p^.) In Fig. 2 we 
see that while with a bare vertex, or even the Ball-Chiu ver- 
sion [7], this function in an Abelian theory is strongly gauge 
dependent, with the Kizilersii-P vertex jst] this is not the case. 
Since in QED the physical coupling is proportional to this 
renormalization function, its independence of the gauge is a 
clear necessity. 
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FIG. 2: Log-log plot of the momentum dependence of the photon 
renormalization function in different covariant gauges defined by ^ = 
-1,0, 1. These are the solutions of the coupled equations shown in 
Fig. 1 with a = g~/4n = 0.2 at p - 10*, for three different ansatze 
for the fermion-boson vertex: (i) bare, (ii) Ball-Chiu [7], (iii) KP 15]. 
Only in the latter case is the result essentially gauge independent 



II. QCD IN THE CONTINUUM 

Studies of the gauge sector of QCD are more complicated. 
It was Baker, Ball and Zachariasen (BBZ) |8] who first used 
the axial gauge to investigate the behaviour of the gluon prop- 
agator They understood the importance of constructing the 
triple gluon vertex to satisfy the relevant Slavnov-Taylor iden- 
tity — the non-Abelian extension of the WGT identity of 
Eq, (3). Axial gauges seem at first sight ideal as they have 
no ghosts and only transverse gluons. However, one never 
gets anything for nothing. In such gauges, the gluon propaga- 
tor depends on two independent functions: coefficients of two 
independent transverse tensors that can be constructed from 
g'^'', the gluon momentum p^' and the axial vector At low- 
est order in perturbation theory, one of these functions is unity, 
and the other is zero. BBZ |8] assumed that this latter func- 
tion remained zero even non-perturbatively, and deduced that 
the gluon propagator behaved like l/p'* at infrared momenta. 
With the potential between infinitely heavy quarks controlled 



by one gluon exchange, this generates a linearly rising "con- 
fining" vector potential. Unfortunately, West |9] showed that 
in axial gauges the gluon could be no more singular than 1 / p^ 
and both gluon functions control the infrared behaviour At- 
tempts to solve the coupled system for these have so far failed 
largely because of the difficulty of dealing consistently with 
the n ■ p singularities non-perturbatively. 

Consequently, attention turned to covariant gauges, in par- 
ticular that of Landau, where of course one has to deal with 
ghosts. In early numerical studies, for instance that by Nick 
Brown and myself ifioll. the ghosts were treated perturb atively, 
just to ensure the ultraviolet behaviour of the gluons was cor- 
rect. However, it was von Smekal and Alkofer [11], who 
showed that ghosts were essential in the infrared too. Indeed, 
they found that the ghosts were singular in the infrared, while 
the gluons were finite or vanishing. They showed these to be 
closely correlated in what is called the scaling solution, with 
the ghost behaving like p^^"^^, while the gluon propagator is 
7?'**"^ as the momentum p ^ Q, with k ^ 0.6. The renormal- 
ization (or dressing) functions (which recall are the propaga- 
tors multiplied by p^) are sketched in Fig. 3. 

For much of hadronic physics the exact behaviour of the 
gluon and ghost propagators in the deep infrared is not rele- 
vant, so we postpone discussion of this till later and instead 
turn to their effect on quarks and then on the hadrons they 
build. In the world of light hadrons, the up and down quarks 
have current masses that are very much less than the scale of 
Aqcd- Even in Abelian theories dynamical mass generation 
can readily occur, provided the interactions are strong enough, 
typically a ^ 7r/3 - 0(1). While perturbatively the vac- 
uum is almost empty, strong long-range interactions change 
its nature. Then particles that appear in the Lagrangian with 
no mass propagate through this medium like heavyweights. 
Since chiral symmetry breaking is such a key feature of the 
low energy strong interaction, modelling this has been a much 
studied aspect of QCD ]fT2lfT3|] . The Schwinger-Dyson equa- 
tion (SDE) for the fermion propagator is shown in Fig. L As 
already remarked, this depends on two functions, the mass 
function M(p^) and the wavefunction renormalization F(p^), 
such that the full propagator is 



Sf(p) 



F(p^) 



P-M{p^)) ' 
where for the bare quantities F - 1 and M - niQ. 



(5) 
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FIG. 3: Gluon and ghost renormalization functions for the scaling 
solution as found by flTll . 
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FIG. 4: Log-log plot of the momentum dependence of the mass func- 
tion M(p^) of Eq. (5) for different values of the current mass of the 
quark that label the lines. The chiral limit of zero current mass is also 
shown. These results from Ref. fl4ll are essentially the same as those 
obtained earlier by Maris and Roberts [ITsIl covering a bigger range 
of masses required for 5 flavours, relevant to the discussion of Fig. 5. 



In QCD, one can study these quark functions alone by mod- 
elling the product of the dressed gluon propagator and the 
quark-gluon vertex (lower graph in Fig. 1), as shown by Maris 
and Tandy 1 12]. Regardless of whether the gluon is enhanced 
in the infrared or suppressed, the key region for physics is the 
momentum region of p ^ A.qcd (with calculations usually 
performed in a momentum subtraction scheme). There the ef- 
fective strength is enhanced, the coupling becomes larger than 
1 and even if mo - a mass is generated. Indeed, one finds 
that for current masses for the up and down quarks defined in 
the perturbative regime to be just a few MeV, a mass of 300- 
400 MeV in the low momentum region is naturally generated, 
as seen in Fig. 4 |14, 15]. Since the SDEs in the continuum 
can be studied for any value of the current mass, one can read- 
ily increase these to 30-100 MeV and so compare with lattice 
results |fl6l] . As shown in iIitIi . these agree remarkably well. 
One of the benefits of working in the continuum is that one 
can also set the quark mass to zero. Then the behaviour of 
the quark propagators differs just a tiny bit from that with 3-5 
MeV (only on a log-log plot like Fig. 4 is this difference to be 
seen). However, in the massless case one can apply the oper- 
ator product expansion to the large momentum behaviour of 
the quark mass function of Fig. 4 and learn that is controlled 
by a qq condensate of -(250 MeV)"*, where the scale is essen- 
tially set by AgcD- Thus it is strong long range correlations 
between quarks and antiquarks that largely breaks chiral sym- 
metry. Moreover, this size of qq condensate agrees with the 
value determined by precision studies of low energy nn final 
state interactions from BNL-E865 [18] and the recent CERN- 
NA48 1 19], and through the Gell-Mann-Oakes-Renner rela- 
tion Il20ll contributes > 90% to the physical pion mass 12 ill . 

This light quark propagator interpolating from small ultra- 
violet current masses to constituent masses below 500 MeV is 



an essential building block for hadronic bound states. Cou- 
pling the Schwinger-Dyson equations with those of Bethe- 
Salpeter allows both the static and dynamical properties of 
light flavoured hadrons to be calculated. This provides a 
framework for studying the relation between hadronic quan- 
tities, for instance, the masses of the lightest pseudoscalar 
and the lightest vector meson. Such studies show the p-mass 
varying with the square of the pion mass, as anticipated from 
the Goldstone nature of the pion. The results naturally fol- 
low those of the lattice for unphysically heavy pions l^22[ |23ll . 
while allowing a natural bridge to the real pion mass, = 
0.02 GeV". Considerable progress has been made with other 
quantum numbers too. These are in general more sensitive to 
details of the approximations beyond the rainbow-ladder for 
the quark-quark scattering kernel, e.g [22*, '24', ^23\. While 
all the calculations discussed here are in Euclidean space, 
or equivalently spacelike momenta in Minkowski space, at- 
tempts to continue to timelike momenta have also been ini- 
tiated |26l]. Similarly detailed dynamical calculations of the 
electromagnetic form-factors of the pion and nucleon have 
been made [27]. They are becoming more realistic and hold 
out the prospect that the momentum behaviour of the up and 
down mass functions shown in Fig. 4 may be amenable to 
experimental test with precision measurements of these form- 
factors |25], planned for the 12 GeV upgrade at JLab. 

Before rushing off and computing more complex hadronic 
effects from these basic quark elements, it is important to re- 
turn to the behavour of the gluons and ghosts themselves that 
are essential for solving the quark SDE of Fig. 1. In fact, this 
was the subject of an exciting and excitable parallel session 
at this workshop. As already mentioned, in the pioneering 
Alkofer, von Smekal, Fischer et al. treatment in the Landau 
gauge ifm l28i l29ll . the gluon is suppressed in the infrared, 
while the ghost is singular. Fig. 3. Importantly, such low 
momentum behaviour is supported by studies [30] of Renor- 
malization Group flow. A singular ghost propagator means it 
dominates almost everywhere it can appear, as discussed by 
Schwenzer et al. [31]. According to [32], the infrared domi- 
nance of the ghost is the key to confinement. At its simplest, 
confinement can be studied non-relativistically by consider- 
ing infinitely heavy quarks. Then the inter-quark potential is 
dominated by one (dressed) gluon exchange with bare ver- 
tices. If this were the case then an infrared vanishing gluon 
would not produce a rising potential at large distances. The 
vertices are expected to be bare because any corrections intro- 
duce additional quark lines and if these have infinite mass they 
are suppressed. The successes of Heavy Quark Effective The- 
ory (HQET) rely on this property. It is here that the singular 
ghost feeds in l32il . The vertex correction with a ghost loop 
brings an infrared singular enhancement to the bare quark- 
gluon vertex, sufficient to generate a confining potential. This 
however comes at a price, the mass function for quarks like 
the b have to behave quite differently from that for the up and 
down quarks. While the b quark is 5 GeV at short distances, 
they find (as shown in the lower part of Fig. 5) its mass func- 
tion is less than 1 GeV at 1 fm. This disagrees with the notions 
that underlie HQET. 
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One may question (y^ whether this behaviour seen in the 
lower graph of Fig. 5 is just an artefact of the approximations 
and truncations made. Let us follow the approach in ll28ll 
and examine the coupled gluon and ghost propagator equa- 
tions, where the interactions, the vertices, are assumed to be 
proportional to their bare structures. As discussed above in 
the context of the photon in QED, Eq. (1), the vector boson 
propagator depends on two functions A and B that gauge 
invariance requires to be equal. Fischer et al. ll28l |29 1 solve 
the equations for A{p^) and this yields the scaling behaviour 
mentioned earlier with k ^ 0.6 . However, explicit calculation 
shows (Fig. 6) this is not the same as B{p^). Though the differ- 
ence is less than a few percent at momenta below 1 GeV ll34ll . 
this makes a substantial difference to the existence of the so- 
lutions. Indeed, demanding these functions are equal, and the 
output gluon is transverse, there is no consistent value for k 
with the approximations used in ll28ll as noted there. 

The key property of the gauge boson propagator is that it 
should be transverse in the Landau gauge both on the right 
and left hand sides of the gluon equation, as well as being 
multiplicatively renormalizable. This means it must not be 
more than logarithmically divergent in the ultraviolet. A non- 
perturbative truncation need not automatically fulfil this. The 
way to regulate the inverse propagator in the Landau gauge for 
a truncated system of Schwinger-Dyson equations is defined 
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(6) 



That this is an appropriate regularisation can be checked by 
considering the contribution of every graph at any order in 
perturbation theory. With dimensional regularisation, every 
graph will be divergent and non-transverse, but when all con- 
tributions at a given order are added they become purely 
logarithmically divergent in the ultraviolet and automatically 
transverse. The correct contribution of every graph to the to- 
tal is ensured by the regularisation given in Eq. (6). This 
is equivalent to determining the gluon's behaviour from the 
function . As noted by Fischer |[2H with the simple in- 

teraction terms used here, a consistent scaling solution then 
requires /c = 1 . However, such a solution, though valid in 
the infrared, has no connection in the ultraviolet to pertur- 
bation theory. Motivated by the Renormalization Group ar- 
guments mentioned earlier [30], Fischer et al. [35] have re- 
cently, within the SDE approach, constructed a more compli- 
cated dressing of the vertices, introducing factors that allow 
the infrared and ultraviolet behaviours to be less tightly linked 
through the shaded region of Fig. 5. Then an infrared scaling 
solution for both A and B functions of Eq. (1) can be formed. 
Nevertheless with a truncated SDE a consistently transverse 
gluon at all momenta requires a regulator like that of Eq. (6). 
Moreover, as illustrated above even for QED, this essential 
regulator does not ensure the boson propagator is physical — 
recall Fig. 2. That requires appropriate interaction terms. In- 
deed, Bloch [36] has emphasised the key role played by the 
two loop graphs in the gluon SDE, which we ignore here, in 
this regard. 

Meanwhile, the Orsay group of Pene et al. Ii37il have long 
investigated the ghost propagator equation in some detail, and 
in particular, the structure of the ghost-gluon vertex for which 
Taylor proved a non-renormalization theorem [38i .39]. They 
have shown how with a vanishing gluon (where k - 0.5 or 
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FIG. 5: Log-log plots of (a) the gluon and ghost renoimalization 
functions, and (b) the quark mass functions from the studies in 
Refs. [28, 32]. The shaded band marks the momentum region that 
largely controls hadron physics. 



FIG. 6: Log-log plot of the momentum dependence of the gluon and 
ghost renormalization functions from an analysis [34] like that of 
\2§]. In the case of the gluon two versions are plotted: short dash 
from p*/A(p^) and the long dash from p'*/B(p^) of Eq. (1). 
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FIG. 7: Linear-log plot of the momentum dependence of a self- 
consistent solution for the gluon and ghost renormalization functions 
with a finite value at = for the ghost function, as described in 
the text m. 

a little bigger) the ghost renormalisation function tends to 
a constant in the infrared, as presented at this meeting by 
Rodriguez-Quintero [%0^. Calculations with my collaborator 
David Wilson |34], confirm the robustness of this conclusion. 
A "scaling" solution is not an automatic consequence of the 
Schwinger-Dyson approach, with its necessary truncations. 

An infrared suppressed gluon more naturally leads to a 
ghost dressing function that is finite as p 0. But can 
such behaviour be found by consistently solving the gluon 
and ghost propagator equations simultaneously? The result 
is sensitive to the exact nature of the ghost-gluon vertex. If 
the ghost loop does dominate the equation for the inverse 
gluon propagator, as in the scaling solution, then Lerche and 
von Smekal |41] have suggested the ghost-gluon vertex itself 
might be transverse to ensure the gluon is consistently trans- 
verse. Of course, a transverse ghost-gluon vertex is too strong 
a requirement, difficult to reconcile with multi-loop orders in 
perturbation theory. The QED experience indicates that a sim- 
pler constraint is required. If one regularises the truncation of 
the gluon propagator equation using Eq. (6), a self-consistent 
outcome for the dressing functions is shown in Fig. 7. The 
result is a "massive" gluon solution, long advocated by Pa- 
pavassiliou and collaborators [42], way back to Cornwall and 
Tiktopoulos ll43ll . The solutions in Fig. 7 have been defined 



by momentum subtraction aX p - Mz, where the ghost and 
gluon renormalization functions are set equal to one. As dis- 
cussed in ll34l[37ll40ll . such a subtraction precludes an infinite 
ghost at zero momentum. The coupling used is equivalent to 
o^siMz) - 0.118 in the MS scheme. For such couplings, a 
massivegluon solution is thus most likely. However Fischer 
et al. 1 35] claim that a scaling solution, as favoured by Renor- 
malization Group flow arguments |30J in the infrared, is still 
possible. So far implementing this within the SDE approach 
introduces ill-constrained vertex dressings fi35il . indicating yet 
more work is required. 

Of course, one must remember that physics of the hadron 
world does not depend on the way ghosts and gluons propa- 
gate over distances of atomic scales, but only those of the size 
of a nucleus — see Fig. 5 for the distance scale. Then the 
difference between a scaling and a massive gluon solution is 
far less marked. Nevertheless the dressed quark propagators, 
using the "massive" gluon of Fig. 7 with a non-singular in- 
teraction in the quark-gluon vertex of Fig. 1, do become like 
those expected of HQET Moreover the ghost and gluon renor- 
malization functions of Fig. 7 are in close accord with a multi- 
plicity of recent lattice results [44, 45 ]. We leave for others (in 
the parallel sessions) to report on the finite volume and lattice 
spacing limitations of such studies. 

What is more the Green's functions of Figs. 7, 4 can readily 
be input into bound state equations to predict the electromag- 
netic form-factors and electro-couplings of excited N*'s that 
are of such physical interest ll25[|27l]. Both the SDE-BSE ap- 
proach in the continuum and lattice studies connect quarks, 
gluons and ghosts to the static and dynamical properties of 
hadrons. In the lattice approach this connection is a "black 
box". However, with the SDE-BSEs it is far less opaque, and 
we can see the inter-related mechanisms at work. It is this 
that makes this a worthwhile long term study. Considerable 
progress has already been made. More is to come. 
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